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Abstract The ground state of dense up and down quark matter under local and global 
charge neutrality conditions with /3-equilibrium has at least four possibilities: 
normal, regular 2SC, gapless 2SC phases, and mixed phase composed of 2SC 
phase and normal components. The discussion is focused on the unusual prop- 
erties of gapless 2SC phase at zero as well as at finite temperature. 
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1. Introduction 

Currently, our knowledge of sufficiently cold and dense matter is very lim- 
ited. There are neither experimental data nor lattice data in this region. From 
the BCS theorem, it is speculated that if the matter is dense enough, the ground 
state of the deconfined quark matter at low temperature will be a color super- 
conductor [1]. Recent studies show that dense quark matter has a rich phase 
structure, see Ref. [2] for reviews. At asymptotically high bayron densites, this 
phenomenon can be studied from first principles [3]. If all three colors of the 
three light quarks participate in the Cooper pairing, the ground state will be in 
the color-flavor-locking (CFL) phase [4]. 

In reality, we are more interested in the intermediate density region, where 
the color superconducting phase may exist in the interior of neutron stars or 
may be created in heavy ion collisions. Unfortunately, we have little knowl- 
edge about this region: we are not sure how the deconfinement and the chiral 
restoration phase transitions happen, how the QCD coupling constant evolves 
and how the strange quark behaves in dense matter, etc. Primarily, our current 
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understanding of the QCD phase structure in this region is based on assump- 
tions. 

In the framework of the bag model [5], or, in general, under the assumption 
that the strange quark mass is small [6], one can exclude the 2SC phase in the 
interior of compact stars when charge neutrality is considered. 

However, there is another possibility, if the strange quark becomes light at 
a larger chemical potential than the u, d quarks, there will be a density re- 
gion where only u, d quarks exist. Here, we focus on the dense quark mat- 
ter composed of only u and d quarks, assuming that the strange quark is too 
heavy to involve in the system. If bulk u, d quark matter exists in the inte- 
rior of compact stars, it should be neutral with respect to electrical as well as 
color charges [6-13]. Also, such matter should remain in /3-equilibrium, i.e., 
the chemical potential for each flavor and color should satisfy the relationship 
Vij,a(3 = (fJ>8ij - HeQij)5 a p + 2/y/3fj,sSij(Ts) a p , where Q and T 8 are genera- 
tors of U(l) em of electromagnetism and the U(l)s subgroup of the color gauge 
group. Satisfying these requirements imposes nontrivial relations between the 
chemical potentials of different quarks. We will see that these requiements 
play very important role in determining the ground state of dense u, d quark 
matter. 

The charge neutrality condition can be satisfied locally [6-11] or globally 
[12, 13]. In the following, we will firstly discuss the homogeneous phase when 
the charge neutrality is satisfied locally, then discuss the mixed phase when the 
charge neutrality condition is satisfied globally. 

2. Local charge neutrality: homogeneous phase 
2.1 Correct way to find the neutral ground state 

To neutralize the electical charge in the homogeneous dense u, d quark 
matter, roughly speaking, twice as many d quarks as u quarks are needed, 
i.e., rid — 2n u , where n u ^ are the number densities for u and d quarks. 
This induces a mismatch between the Fermi surfaces of pairing quarks, i.e., 
fJ-d — = fJ-e = 2<5/x, where fi e is the electron chemical potential. 

To get the ground state of the system, we need to know the thermodynamical 
potential. For simplicity, we use Nambu-Jona-Lasinio (NJL) model [14] to 
describe 2-flavor quark matter, 

C = qij^d^q + G s {qqf + {qi^rqf 

+G D [(iq c ee b l5 q)(iqee b 75 q c )] , (1) 

where q c = Cq T is the charge-conjugate spinor and C = ij 2 j is the charge 
conjugation matrix. The quark field q = qi a is a four-component Dirac spinor 
that carries flavor (i = 1, 2) and color (a = 1,2,3) indices. The Pauli matrices 
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are denoted by f = (r 1 ^ 2 ^ 3 ), while {e) tk = e lk and (e b ) afS = e afSb are 
the antisymmetric tensors in the flavor and color spaces, respectively. We also 
introduce a momentum cutoff A, and two independent coupling constants in 
the scalar quark- antiquark and scalar diquark channels, Gs and Gd- We define 
77 = Gd/Gs, and 77 = 0.75 is the value from the Fierz transformation, but in 
principle, 77 is a free parameter. 

In the mean-field approximation, the thermodynamical potential for u, d 
quarks in /3-equilibrium with electrons takes the form [7, 9]: 

= !i °-i^("« +2 ^"- + ^ T4 ) + ^ 

+ ^-?/|^+ 2 ™(^- E " /T )]> <» 

where fi is a constant added to make the pressure of the vacuum zero, and the 
electron mass was taken to be zero, which is sufficient for the purposes of the 
current study. The sum in the second line of Eq. (2) runs over all (6 quark and 6 
antiquark) quasi-particles. The explicit dispersion relations and the degeneracy 
factors of the quasi-particles read 

E± = E(p)±fi ub , [xl] (3) 

E± = E(p)±fi db , [xl] (4) 

E± ± = E±(p)±6fi. [x2] (5) 

Here we introduced the following shorthand notation: E(p) = ^p 2 + m 2 and 
E±(p) = y/[E(p) ± fi\ 2 + A 2 with p, = fi — fi e /6 + n 8 /3. 

If a macroscopic chunk of quark matter is created in heavy ion collisions or 
exists inside the compact stars, it must be in color singlet. So in the following 
discussions, color charge neutrality condition is always satisfied. 

Now, we discuss the role of electrical charge neutrality condition. If a 
macroscopic chunk of quark matter has nonzero net electrical charge density 
tiq, the total thermodynamical potential for the system should be given by 

= ^Coulomb + ^u,d,e; (6) 

where Slcoulomb ~ n Q^ 2 ^ 3 (V i s the volume of the system) is induced by 
the repulsive Coulomb interaction. The energy density grows with increasing 
the volume of the system, as a result, it is almost impossible for matter inside 
stars to remain charged over macroscopic distances. So the bulk quark matter 
should also satisfy electrical neutrality condition, thus Qcouiomb\n Q =o = 0, 
and £l u ,d,e\n Q =o is on the neutrality line. Under the charge neutrality condition, 
the total thermodynamical potential of the system is Q\ nQ=0 = VL u ^ e \n Q =a- 
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Figure 1. The effective potential as a function of the diquark gap A calculated at a fixed value 
of the electrical chemical potential /i e = 148 MeV (dashed line), and the effective potential 
defined along the neutrality line (solid line). The results are plotted for fx = 400 MeV with 
t) = 0.75. 



Here, we want to emphasize that: The correct way to find the ground state of 
the homogeneous neutral u, d quark matter is to minimize the thermodynam- 
ical potential along the neutrality line 0| nQ= o = ^u,d,e\n Q =o, not like in the 
flavor asymmetric quark system, where (3-equilibrium is required but fj, e is a 
free parameter, and the ground state is determined by minimizing the thermo- 
dynamical potential fl u ,d,e- 

From Figure 1, we can see the difference in determining the ground state 
for a charge neutral system and for a flavor asymmetric system. In Figure 1, 
at a given chemical potential \i = 400 MeV and 77 = 0.75, the thermody- 
namical potential along the charge neutrality line Q\ nQ= o as a function of the 
diquark gap A is shown by the solid line. The minimum gives the ground state 
of the neutral system, and the corresponding values of the chemical potential 
and the diquark gap are [i e = 148 MeV and A = 68 MeV, respectively. If 
we switch off the charge neutrality conditions, and consider the flavor asym- 
metric u, d quark matter in 0— equilibrium [15], the electrical chemical poten- 
tial n e becomes a free parameter. At a fixed [i e = 148 MeV and with color 
charge neutrality, the thermodynamical potential is shown as a function of the 
diquark gap by the dashed line in Figure 1. The minimum gives the ground 
state of the flavor asymmetric system, and the corresponding diquark gap is 
A = 104 MeV, but this state has positive electrical charge density, and cannot 
exist in the interior of compact stars. 



Neutral Dense Quark Matter 



5 



A (MeV) 




Figure 2. The graphical representation of the solution to the charge neutrality conditions 
(thick dash-dotted line) and the solution to the gap equation for three different values of the 
diquark coupling constant (thick solid and dashed lines). The intersection points represent the 
solutions to both. The thin solid line divides two qualitatively different regions, A < <5/i and 
A > S[i. The results are plotted for /i = 400 MeV and three values of diquark coupling 
constant Gd = r/Gs with rj = 0.5, r\ — 0.75, and r\ = 1.0. 

2.2 t) dependent solutions 

In the last subsection, by looking for the minimum of the thermodynamic 
potential along the charge neutrality line, we found the ground state for the 
charge neutral u, d quark system. 

Equivalently, the neutral ground state can also be determined by solving 
the diquark gap equation together with the charge neutrality conditions. We 
visualize this method in Figure 2, with color neutrality always satisfied, at a 
given chemical potential fi = 400 MeV. The nontrivial solutions to the diquark 
gap equation as functions of the electrical chemical potential \i e are shown by 
a thick-solid line {r\ = 0.75), a long-dashed line (r? = 1.0), and a short-dashed 
line (77 = 0.5). It is found that for each 77, the solution is divided into two 
branches by the thin-solid line A = 5fi, and the solution is very sensitive to 
77. Also, there is always a trivial solution to the diquark gap equation, i.e., 
A = 0. The solution of the charge neutrality conditions is shown by a thick 
dash-dotted line, which is also divided into two branches by the thin-solid line 
A = Sfi, but the solution of the charge neutrality is independent of rj. 

The cross-point of the solutions to the charge neutrality conditions and the 
diquark gap gives the solution of the system. We find that the neutral ground 
state is sensitive to the coupling constant Gd = f]Gs in the diquark channel. 
In the case of a very strong coupling (e.g., r\ = 1.0 case), the charge neutrality 
line crosses the upper branch of the solution to the diquark gap, the ground 
state is a charge neutral regular 2SC phase with A > 8/j,. In the case of weak 
coupling (e.g., rj = 0.5), the charge neutrality line crosses the trivial solution 
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Figure 3. The quasi-particle dispersion relations at low energies in the 2SC phase (left panel) 
and in the g2SC phase (right panel). 

of the diquark gap, i.e., the ground state is a charge neutral normal quark matter 
with A = 0. The regime of intermediate coupling (see, e.g., r] = 0.75 case) 
is most interesting, the charge neutrality line crosses the lower branch of the 
solution of the diquark gap. We will see that this phase is a gapless 2SC (g2SC) 
phase with A < 8/j,, which is different from the regular 2SC phase, and has 
some unusual properties. 

2.3 g2SC phase 

In this subsection, we will explain why we call the color superconducting 
phase with A < 5/j, the g2SC phase, and we will show some special properties 
of this phase. 

Quasi-particle spectrum 

It is instructive to start with the excitation spectrum in the case of the ordi- 
nary 2SC phase when S/j, = 0. With the conventional choice of the gap point- 
ing in the anti-blue direction in color space, the blue quarks are not affected by 
the pairing dynamics, and the other four quarsi-particle excitations are linear 
superpositions of u r:9 and d r , g quarks and holes. The quasi-particle is nearly 
identical with a quark at large momenta and with a hole at small momenta. We 
represent the quasi-particle in the form of Q(quark, hole), then the four quasi- 
particles can be represented explicitly as Q(u r , d g ), Q(u g , d r ), Q(d r , u g ) and 
Q(d g ,u r ). When 5/j, = 0, the four quasi-particles are degenerate, and have a 
common gap A. 



Neutral Dense Quark Matter 



7 



If there is a small mismatch (<5/x < A) between the Fermi surfaces of the 
pairing u and d quarks, the excitation spectrum will change. For example, 
we show the excitation spectrum of Q(u r , d g ) and Q(d g , u r ) in the left panel 
of Figure 3. We can see that S/j, induces two different dispersion relations, 
the quasi -particle Q(d g ,u r ) has a smaller gap A — //, and the quasi-particle 
Q(u r , d g ) has a larger gap A+/i. This is similar to the case when the mismatch 
is induced by the mass difference of the pairing quarks [16]. 

If the mismatch 6/j, is larger than the gap parameter A, the lower disper- 
sion relation for the quasi-particle Q(d g ,u r ) will cross the zero-energy axis, 
as shown in the right panel of Figure 3. The energy of the quasi-particle 
Q(d g ,u r ) vanishes at two values of momenta p = pT and p = ^+ where 
/r* 1 = p, ± y/ (<5/Lt) 2 — A 2 . This is why we call this phase gapless 2SC (g2SC) 
phase. An unstable gapless CFL phase has been found in Ref. [16], and a sim- 
ilar stable gapless color superconductivity could also appear in a cold atomic 
gas or in u, s or d, s quark matter when the number densities are kept fixed 
[17]. 

As one would expect, far outside the pairing region, p ~ fi, the quasi- 
particle dispersion relations are similar to those in the 2SC phase. Also, around 
p ~ jl, the quasi-particle Q(u r , d g ) resembles the dispersion relations with that 
in the regular 2SC phase. The most remarkable property of the quasi-particle 
spectra in the g2SC phase is that the low energy excitations (E <C Sfi — A) 
are very similar to those in the normal phase represented by solid lines. The 
only difference is that the values of the chemical potentials of the up and down 
quarks fi ur = p, ug and Hdg = ^dr are replaced by the values ix~ and fi + , 
respectively. This observation suggests, in particular, that the low energy (large 
distance scale) properties of the g2SC phase should look similar to those in the 
normal phase. 

Finite temperature properties 

In a superconducting system, when one increases the temperature at a given 
chemical potential, thermal motion will eventually break up the quark Cooper 
pairs. In the weakly interacting Bardeen-Copper-Schrieffer (BCS) theory, the 
transition between the superconducting and normal phases is usually of second 
order. The ratio of the critical temperature T C BCS to the zero temperature value 
of the gap Aq cs is a universal value [18] 

T BCS e 7 B 

r B CS = ^jp = — ~ °- 567 ' W 

where je ~ 0.577 is the Euler constant. In the conventional 2SC phase of 
quark matter with equal densities of the up and down quarks, the ratio of the 
critical temperature to the zero temperature value of the gap is also the same 
as in the BCS theory [19]. In the spin-0 color flavor locked phase as well as in 
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Figure 4. The temperature dependence of the diquark gap (solid line) and the value of Sfi = 
/x e /2 (dashed line) in neutral quark matter. For comparison, the diquark gap in the model with 
/i e = and [is = is also shown (dash-dotted line). The results are plotted for /i = 400 MeV 
and rj = 0.75. 

the spin-1 color spin locked phase, on the other hand, this ratio is larger than 
the BCS ratio by the factors 2 1 / 3 and 2 2 / 3 , respectively. These deviations are 
related directly to the presence of two different types of quasi-particles with 
nonequal gaps [20]. 

For the g2SC phase, the typical results for the default choice of parameters 
\i = 400 MeV and rj = 0.75 are shown in Figure 4. Both the values of the 
diquark gap (solid line) and the mismatch parameter 5[i = /U e /2 (dashed line) 
are plotted. One very unusual property of the shown temperature dependence 
of the gap is a nonmonotonic behavior. Only at sufficiently high temperatures, 
the gap is a decreasing function. In the low temperature region, T < 10 MeV, 
however, it increases with temperature. For comparison, in the same figure, the 
diquark gap in the model with fi e = and /ig = is also shown (dash-dotted 
line). This latter has the standard BCS shape. 

Another interesting thing regarding the temperature dependences in Figure 4 
appears in the intermediate temperature region, 22.5 < T < 37 MeV. By com- 
paring the values of A(T) and 5[i in this region, we see that the g2SC phase 
is replaced by a "transitional" 2SC phase there. Indeed, the energy spectrum 
of the quasi-particles even at finite temperature is determined by the same re- 
lations in Eqs. (3) and (5) that we used at zero temperature. When A > S/j,, 
the modes determined by Eq. (5) are gapped. Then, according to our standard 
classification, the ground state is the 2SC phase. 

It is fair to say, of course, that the qualitative difference of the g2SC and 2SC 
phases is not so striking at finite temperature as it is at zero temperature. This 
difference is particularly negligible in the region of interest where temperatures 
22.5 < T < 37 MeV are considerably larger than the actual value of the 
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smaller gap, A — 5fi. However, by increasing the value of the coupling constant 
slightly, the transitional 2SC phase can be made much stronger and the window 
of intermediate temperatures can become considerably wider. In either case, 
we find it rather unusual that the g2SC phase of neutral quark matter is replaced 
by a transitional 2SC phase at intermediate temperatures which, is replaced by 
the g2SC phase again at higher temperatures. 

It appears that the temperature dependence of the diquark gap is very sensi- 
tive to the choice of the diquark coupling strength 7] = Gd/Gs in the model 
at hand. This is not surprising because the solution to the gap equation is very 
sensitive to this choice. The resulting interplay of the solution for A with 
the condition of charge neutrality, however, is very interesting. This is demon- 
strated by the plot of the temperature dependence of the diquark gap calculated 
for several values of the diquark coupling constant in Figure 5. 

The most amazing are the results for weak coupling. It appears that the gap 
function could have sizable values at finite temperature even if it is exactly zero 
at zero temperature. This possibility comes about only because of the strong 
influence of the neutrality condition on the ground state preference in quark 
matter. Because of the thermal effects, the positive electrical charge of the 
diquark condensate is easier to accommodate at finite temperature. We should 
mention that somewhat similar results for the temperature dependence of the 
gap were also obtained in Ref. [21] in a study of the asymmetric nuclear matter, 
and in Ref. [22] when number density was fixed. 

The numerical results for the ratio of the critical temperature to the zero 
temperature gap in the g2SC case as a function of the diquark coupling strength 
rj = Gd/Gs are plotted in Figure 6. The dependence is shown for the most 
interesting range of values of rj = Gd/ Gs, 0.68 < rj < 0.81, which allows the 
g2SC stable ground state at zero temperature. When the coupling gets weaker 




in this range, the zero temperature gap vanishes gradually. As we saw from 
Figure 5, however, this does not mean that the critical temperature vanishes 
too. Therefore, the ratio of a finite value of T c to the vanishing value of the gap 
can become arbitrarily large. In fact, it remains strictly infinite for a range of 
couplings. 

3. Global charge neutrality: mixed phase 

We have discussed the homogeneous 2-flavor quark matter when charge 
neutrality conditions are satisfied locally, and found that the local charge neu- 
trality conditions impose very strong constraints on determining the ground 
state of the system. 

On the other hand, one can construct mixed phase when charge neutrality 
conditions are satisfied globally. Inside mixed phases, the charge neutrality is 
satisfied "on average" rather than locally. This means that different compo- 
nents of mixed phases may have non-zero densities of conserved charges, but 
the total charge of all components still vanishes. In this case, one says that the 
local charge neutrality condition is replaced by a global one. There are three 
possible components: (i) normal phase, (ii) 2SC phase, and (iii) g2SC phase. 

The pressure of the main three phases of two-flavor quark matter as a func- 
tion of the baryon and electrical chemical potentials is shown in Figure 7 at 
T) = 0.75. In this figure, we also show the pressure of the neutral normal quark 
and gapless 2SC phases (two dark solid lines). The surface of the g2SC phase 
extends only over a finite range of the values of fi e . It merges with the pres- 
sure surfaces of the normal quark phase (on the left) and with the ordinary 2SC 
phase (on the right). 




Figure 7. At r\ — 0.75, pressure as a function of /i = /j.b/3 and [i e for the normal and color 
superconducting quark phases. The dark solid lines represent two locally neutral phases: (i) the 
neutral normal quark phase on the left, and (ii) the neutral gapless 2SC phase on the right. The 
appearance of the swallowtail structure is related to the first order type of the phase transition in 
quark matter. 



It is interesting to notice that the three pressure surfaces in Figure 7 form 
a characteristic swallowtail structure. As one could see, the appearance of 
this structure is directly related to the fact that the phase transition between 
color superconducting and normal quark matter, which is driven by changing 
parameter fj, e , is of first order. In fact, one should expect the appearance of a 
similar swallowtail structure also in a self-consistent description of the hadron- 
quark phase transition. Such a description, however, is not available yet. 

From Figure 7, one could see that the surfaces of normal and 2SC quark 
phases intersect along a common line. This means that the two phases have the 
same pressure along this line, and therefore could potentially co-exist. More- 
over, as is easy to check, normal quark matter is negatively charged, while 2SC 
quark matter is positively charged on this line. This observation suggests that 
the appearance of the corresponding mixed phase is almost inevitable. 

Let us start by giving a brief introduction into the general method of con- 
structing mixed phases by imposing the Gibbs conditions of equilibrium [23, 
24]. From the physical point of view, the Gibbs conditions enforce the me- 
chanical as well as chemical equilibrium between different components of a 
mixed phase. This is achieved by requiring that the pressure of different com- 
ponents inside the mixed phase are equal, and that the chemical potentials (pt 
and fi e ) are the same across the whole mixed phase. For example, in relation 
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Figure 8. At r\ — 0.75, pressure as a function of /i = and /j, e for the normal and 

color superconducting quark phases (the same as in Figure 7, but from a different viewpoint). 
The dark solid line represents the mixed phase of negatively charged normal quark matter and 
positively charged 2SC matter. 



to the mixed phase of normal and 2SC quark matter, these conditions read 

P^Gu^e) = pCa^W), (8) 
// = ^=^ 2SC \ (9) 

Me = M fQ) =/i (25C) (10) 

It is easy to visualize these conditions by plotting the pressure as a function 
of chemical potentials (/i and fi e ) for both components of the mixed phase. 
This is shown in Figure 8. As should be clear, the above Gibbs conditions 
are automatically satisfied along the intersection line of two pressure surfaces 
(dark solid line in Figure 8). 

Different components of the mixed phase occupy different volumes of space. 
To describe this quantitatively, we introduce the volume fraction of normal 
quark matter as follows: X-2SC = Vnq/V (notation Xb means volume frac- 
tion of phase A in a mixture with phase B). Then, the volume fraction of the 
2SC phase is given by x'n'q = (1 — X2sc)- F rom the definition, it is clear that 

o<xS<i. 

The average electrical charge density of the mixed phase is determined by 
the charge densities of its components taken in the proportion of the corre- 
sponding volume fractions. Thus, 

4 MP) = iSM'^M + (1 - X N 2 s Q c)nf SC) {^^). (ID 
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If the charge densities of the two components have opposite signs, one can im- 
pose the global charge neutrality condition, n{ MP ^ = 0. Otherwise, a neutral 
mixed phase could not exist. In the case of quark matter, the charge density 
of the normal quark phase is negative, while the charge density of the 2SC 
phase is positive along the line of the Gibbs construction (dark solid line in 
Figure 8). Therefore, a neutral mixed phase exists. The volume fractions of its 
components are 



(2SC) 
i-e 

n r ^ - ni 



A2SC - (2SC) (NQ)' ylA> 



^2SC _ -, NQ _ n e 

Xnq = 1 - X 2 sc - ~Jnq] (2scy {U) 

Tie fie 

After the volume fractions have been determined from the condition of the 
global charge neutrality, we could also calculate the energy density of the cor- 
responding mixed phase, 

e {MP) = X^ (iVQ) (^e) + (1 - X&%) £ < MC W). (14) 

This is essentially all that we need in order to construct the equation of state of 
the mixed phase. 

So far, we were neglecting the effects of the Coulomb forces and the surface 
tension between different components of the mixed phase. In a real system, 
however, these are important. In particular, the balance between the Coulomb 
forces and the surface tension determines the size and geometry of different 
components inside the mixed phase. 

In our case, nearly equal volume fractions of the two quark phases are likely 
to form alternating layers (slabs) of matter. The energy cost per unit volume to 

produce such layers scales as <r 2 / 3 (n e 2SC ^ — n e N ^) 2 ^ 3 where a is the surface 
tension [25]. Therefore, the quark mixed phase is a favorable phase of matter 
only if the surface tension is not too large. Our simple estimates show that 
&max < 20 MeV/fm 2 . However, even for slightly larger values, 20 < a < 50 
MeV/fm 2 , the mixed phase is still possible, but its first appearance would occur 
at larger densities, 3po < ps < 5po- The value of the maximum surface 
tension obtained here is comparable to the estimate in the case of the hadronic- 
CFL mixed phase obtained in Ref. [26]. The thickness of the layers scales 

as (j 1 / 3 (ne 2 ' 5C7 ' ) — n e N ^)~ 2 / 3 [25], and its typical value is of order 10 fm in 
the quark mixed phase. This is similar to the estimates in various hadron- 
quark and hadron-hadron mixed phases [25, 26]. While the actual value of the 
surface tension in quark matter is not known, in this study we assume that it is 
not very large. Otherwise, the homogeneous gapless 2SC phase should be the 
most favorable phase of nonstrange quark matter [9]. 
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Under the assumptions that the effect of Coulomb forces and the surface 
tension is small, the mixed phase of normal and 2SC quark matter is the most 
favorable neutral phase of matter in the model at hand with 77 = 0.75. This 
should be clear from observing the pressure surfaces in Figs. 7 and 8. For a 
given value of the baryon chemical potential \i = /ie/3, the mixed phase is 
more favorable than the gapless 2SC phase, while the gapless 2SC phase is 
more favorable than the neutral normal quark phase. 

4. Conclusion 

Dense it, d quark matter under local and global charge neutrality conditions 
in /3-equilibrium has been discussed. 

Under local charge neutrality condition, the homogeneous neutral ground 
state is sensitive to the coupling constant in the diquark channel, it will be in the 
regular 2SC phase when the coupling is strong, in the normal phase when the 
coupling is weak, and in the g2SC phase in the case of intermediate coupling. 
The low energy quasi-particle spectrum in g2SC phase contains four gapless 
and only two gapped modes, and this phase has rather unusual properties at 
zero as well as at finite temperature. 

Under global charge neutrality condition, assuming that the effect of Coulomb 
forces and the surface tension is small, one can construct a mixed phase com- 
posed of positive charged 2SC phase and negative charged normal quark mat- 
ter. 
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